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Abstract 

We study lattice constructions of gapped fermionic phases of mat¬ 
ter. We show that the construction of fermionic Symmetry Protected 
Topological orders by Gu and Wen has a hidden dependence on a 
discrete spin structure on the Euclidean space-time. The spin struc¬ 
ture is needed to resolve ambiguities which are otherwise present. An 
identical ambiguity is shown to arise in the fermionic analog of the 
string-net construction of 2D topological orders. We argue that the 
need for a spin structure is a general feature of lattice models with 
local fermionic degrees of freedom and is a lattice analog of the spin- 
statistics relation. 


1 Introduction and summary 

1.1 Bosonic and fermionic gapped phases 

In condensed matter physics, topological phases of matter are often 
defined as equivalence classes of local gapped bosonic Hamiltonians, 
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usually defined on a lattice, which can be deformed into each other 
without ever becoming gapless mm- The notion of topological phase 
can be enriched by imposing additional constraints on the theories, 
such as a choice of global symmetry preserved by all the Hamiltonians. 

On the other hand, topological quantum field theorie^ can be thought 
of as describing the far infrared behavior of gapped bosonic quantum 
field theories (see e.g. section 4 of [3] or the monograph HD- 

There is a close relation between topological phases of matter and 
topological quantum field theories, which can be thought of as a map 
from a topological phase of matter to the TQFT which encodes the low 
energy continuum limit of the corresponding Hamiltonian. In princi¬ 
ple, one may imagine the map being many-to-one; it is not obvious 
that two local gapped Hamiltonians which map to the same TQFT 
will always be deformable into each other. Still, in practice we do not 
know of any observable which can distinguish two phases of matter, 
but cannot be formulated in terms of the TQFT data. El 

In condensed matter physics, one also encounters the notion of a 
fermionic topological phase of matter, defined as an equivalence class 
of local gapped Hamiltonians which can involve fermionic degrees of 
freedom urn- Perhaps surprisingly, some fermionic phases of matter 
are not expected to admit a purely bosonic realization. This is ex¬ 
pected to be due to the difference in the notion of locality for bosonic 
and fermionic systems. Intuitively, if we partition a bosonic system in 
two parts, the total Hilbert space factors uniquely in the tensor prod¬ 
uct of the Hilbert spaces for the two parts. If we partition a fermionic 
system, though, the factorization has an intrinsic ambiguity, as ob¬ 
servables in the tensor product of the Hilbert spaces for the two parts 
are defined up to a sign in the sector where both factors have odd 
fermion number. 

^This is a somewhat looser notion of TQFT compared to some formal definitions. For 
example, we consider Chern-Simons theory to be a TQFT, even though it has a partition 
function which depends on a choice of metric on space-time. In other words, we allow the 
stress tensor to be non-zero, but proportional to the identity operator. 

^It is also conceivable, perhaps, that some topological phase of matter may not give rise 
to a TQFT at low energy, i.e. that some anomaly/obstruction may prevent the definition 
of TQFT amplitudes on general manifolds in terms of the Hamiltonian data. But in 
all cases known to us one circumvent such obstructions by postulating that the TQFT 
depends on some additional geometric data, such as metric or framing. 
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1.2 Spin structure dependence 

In unitary quantum field theory, fermions are naturally spinors and 
thus the low energy physics of a gapped fermionic theory is a spin- 
TQFT: a topological field theory defined on manifolds which can be 
equipped with a spin structure, whose correlation functions possibly 
depend on the choice of spin structure. The purpose of this paper is 
to explore the relation between fermionic topological phases of matter 
and spin-TQFTs. It is not obvious that such a relation should exist, as 
a lattice Hamiltonian involving fermionic degrees of freedom is usually 
written down without any reference to a spin structure on the manifold 
which is discretized by the lattice. One also cannot appeal to the 
spin-statistics relation, because the lattice destroys Lorenz and even 
rotational invariance which are the conditions of the spin-statistics 
theorem, d 

The first step of our analysis is to look carefully at the fermionic 
SPT phases constructed by Gu and Wen in [7]. We find that the pre¬ 
scription used to define the partition function of such theories runs into 
an obstruction if applied to space-time manifolds of general topology, 
unless the second Stiefel-Whitney class [^ 2 ] of the manifolds vanishes, 
i.e. the manifold admits a spin structure. If the manifolds admits a 
spin structure, the obstruction can be eliminated, but the final answer 
will depend on the choice of spin structure t]. In other words, these 
fermionic SPT phases define (invertible) spin-TQFTs. 

Next, we look at other known constructions of fermionic phases 
of matter which are expected to admit a state-sum-like definition of 
their partition function: the construction of fermionic toric code in 
[8] and the general fermionic Turaev-Viro construction in [9]. These 
references focus on the construction of a fixed-point Hamiltonian and 
wave-function for these fermionic phases of matter, rather than a par¬ 
tition function. It is straightforward, though, to assemble the same 
ingredients into a partition sum, borrowing some ideas from the Gu- 
Wen fermionic SPT phase construction. Again, we find an obstruction 
to define the partition sum unless the space-time manifolds admits a 

^Taking the continuum limit and then applying the spin-statistics relation does not 
ameliorate the problem. The TQFT itself is, of course, Lorentz invariant, but the spin- 
statistics relation is a property of Lorenz invariant particle excitations. The continuum 
limit from the lattice theory to the low-energy TQFT only concerns the ground states 
of the system. A priori, massive excitations above these ground states do not need to 
transform properly under the Lorentz group. 
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spin structure, in which case one can remove the obstruction and de¬ 
fine a well defined partition function which depends on the choice of 
spin structure r/. Thus these fermionic phases of matter are associated 
to spin-TQFTs. 

We can describe the obstruction schematically here, referring the 
reader to sections [2] and El for further details. State-sum models as¬ 
semble the partition function from a triangulation of the space-time 
manifold X: each simplex is associated to some tensor in the tensor 
product of vector spaces associated to the faces and the legs of these 
tensors are contracted together as the simplices are glued along the 
corresponding faces of the triangulation. In a fermionic model, the 
vector spaces may be Grassmann-odd and Koszul signs occur when 
re-organizing and contracting the factors of the tensor products. 

These Koszul signs, arising from the anti-commutation of fermionic 
variables, are of course a key element of the problem. The non-local 
nature of these signs is precisely what should allow these fermionic 
phases of matter to be distinct from any bosonic phase. In order for the 
partition sum to be invariant under local changes in the triangulation 
of the manifold, one needs to cancel the change in the Koszul signs 
agains the change in the local data attached to the simplices. The 
obstruction arises precisely when this cancellation is not possible. 

We can express the obstruction neatly by encoding the fermion 
number of the vector spaces attached to faces in a Z 2 -valued {d — 1)- 
cochain f3d-i- The cochain /3d-i is actually a cocycle, as the total 
fermion number of the tensors attached to simplices is even. It is useful 
to decompose the partition sum into a sum of terms Z[X, Pd-i], which 
contain the parts of the state sum due to states of fermion number 
Pd-l- 

We can encode a general change of triangulation of X into a tri¬ 
angulation of the {d + l)-dimensional manifold X x [0,1]. Intuitively, 
we are gluing a sequence of {d -|- l)-dimensional simplices on top of 
our initial triangulation to get the final triangulation. We find that 
the triangulation invariance of the partition function is obstructed by 
some irreducible sign mismatch, which can be written schematically 
as 

Here W 2 is a 2-cocycle with values in Z 2 representing the second Stiefel- 
Whitney class of X x [0,1] and Pd-i is a lift to X x [0,1] of the cocycle 
fdd-i- If the cohomology class of \w‘ 2 \ is non-trivial and the theory in¬ 
volves choices of fermion numbers f^d-i which are non-trivial in coho- 
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mology, this sign mismatch cannot be absorbed by a redefinition of the 
local part of the partition function. This prevents us from constructing 
a well-defined partition sum and ruins the state-sum construction. 

If we restrict X to be a spin manifold then W 2 {X) is exact and we 
can write W 2 = Sr] for some 1-cochain t], which represents a choice of 
spin structure. This allows us to thus cancel the obstruction [T] by the 
variation of a local term 

(_l)/x^U/3d-i (2) 

so that the improved state sum 

Z[X,r,] = (3) 

Pd-l 

is fully invariant under changes of triangulations and defines a good 
theory. This theory is a spin-TQFT: it can only be defined on a spin 
manifold and depends on a choice of spin structure. 

In sections [3] and m we will look in further detail at the properties 
of the Koszul signs which occur in the state sum. The definition the 
partition function requires specific choices of how to order the factors 
in the tensor product associated to each simplex, and the two factors 
in the contraction of vector spaces at each face. Given some ordering 
choices 11, the permutations of the vector spaces involved in the state 
sum will produce some overall Koszul sign aYi{X, /3rf_i), which depends 
only on the triangulation, on 11 and on (5d-i- 

The choice of order 11 can be given independently of the other data 
in the state sum. The combined sign 

zn[X,r,,/3rf_i] =^7n(X,/3rf_i)(-l)/v^u/3,_i ( 4 ) 

appears to be a very useful object, which captures the intrinsically 
fermionic part of the full partition function. From now on we will 
drop the subscript 11. Our formulae will refer to the specific choice of 
order used in the Gu-Wen definition of fermionic SPT phases [7]. We 
will comment briefly on other choices of order in section [3l 

We can think about z[X,ri, as defining an effective action 

for a {d — I)-form Z 2 gauge field with a very specific anomaly, or a 
very simple invertible spin-TQFT equipped with an anomalous 
{d — 2)-form Z 2 global symmetry. 

Under changes of triangulation, z[X,r], /3d-i] changes by another 
interesting cocycle, the Steenrod square of /3d-i- 
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We refer to appendix [B] for the explicit definition of the higher cup 
products Uq. Their basic property is 

AUaB + BUaA = 5{A Ua+i B) + 5A Ua+i B + A Ua+i 5B ( 6 ) 

with Uq = U. 

Under gauge transformations, we find the precise form of the ’t 
Hooft anomaly 

Z[X, r?, Pd-l+SXd- 2 ] = z[X, rj, /3d_i](-l)/x /3Ud-3A+AU,_3/3+AU,_35A+AU,_4A 

(7) 

Although z[X, r], Pd-i] does not appear to admit a d-dimensional 
bosonic description, we also find that it is a quadratic refinement of a 
bosonic pairing: 

( 8 ) 

Finally, if X is a boundary of a compact oriented {d + l)-manifold Y 
and d > 2, we find an explicit WZW-like expression for z[X,r], /3d-i]- 

z[X,r],Pd-l] = (-l)/x^U/3,_i+/ySg2[/3d_i]+«;2U/3d_i 

Here we use the fact that for d > 2 the cocycle Pd-i can be extended to 
Y. The action is independent of the choice of Y or of the way fid-i is 
extended from X too Y because the expression Sq'^[Pd-i] +rc 2 U Pd-i 
is exact for closed oriented Y . This formula is particularly useful 
for d = 3, since any closed oriented 3-manifold X is a boundary of a 
compact oriented 4-manifold Y. 

With a bit of extra work, we can rewrite the partition function 
Z[X, r]] of our spin-TQFT as the partition function of a (d — l)-form 
Z 2 gauge theory 

Z[X, v] = Yl Pd-i]z[X, rj, Pd-i] (10) 

Pd-i 


where the gauge fields are coupled to two two sets of degrees of free¬ 
dom: a standard bosonic TQFT equipped with a (d—2)-form Z 2 global 
symmetry and partition function Z[X, I3d-i], and the spin-TQFT Kd- 
The bosonic theory associated to Z[X,f3d-i] must have a ’t Hooft 
anomaly which cancels the ’t Hooft anomaly of Kd, controlled by 
Sq^lPd-i]- 

In order to make contact with concepts which are more famil¬ 
iar in condensed matter physics, it is useful to replace the notion of 
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a TQFT with an anomalous global symmetry with the notion of a 
gapped boundary condition for a (d+ l)-dimensional SPT phase, pro¬ 
tected by a (d — 2)-form Z 2 global symmetry, with partition function 

Then Z\X, Pd-i] defines a bosonic gapped boundary condition for the 
{d+ l)-dimensional SPT phase, while z[X, rj, fid-i] defines a fermionic 
gapped boundary condition. The original spin TQFT can be recovered 
by gauging the {d — 2)-form Z 2 global symmetry on a slab, with one 
of these boundary conditions at either end. 

In section [7] we will argue that this construction has a close rela¬ 
tion to the notion of fermionic anyon condensation. The generators 
of a non-anomalous {d — 2)-form global symmetry G are loop observ¬ 
ables which can be thought as worldlines of bosonic quasi-particles 
which fuse accordingly to the group law of G. Gauging the (d — 2)- 
form symmetry is equivalent to proliferating these quasi-particles in 
correlation functions and should correspond to the standard notion of 
anyon condensation, at least in three space-time dimensions. 

We will argue that the generators of a (d — 2)-form Z 2 global sym¬ 
metry with the t’Hooft anomaly described above, instead, behave as 
fermionic quasi-particles. The challenge to define a fermionic analogue 
of the standard anyon condensation is mapped into the problem of 
gauging such anomalous symmetry. The kernel spin-TQFT Kd offers 
a solution to the problem: given some TQFT with fermionic quasi¬ 
particles we want to condense, we can tensor it with Kd to cancel the 
anomaly and gauge the (d — 2)-form Z 2 global symmetry. Essentially, 
the kernel spin-TQFT Kd uses the spin structure information to pro¬ 
vide some extra signs which make the “fermionic” anyon condensation 
meaningful. This approach to fermionic anyon condensation appears 
to be closely related to work in progress by K. Walker [25j . 

State-sum constructions of spin-TQFTs in two dimensions have 
been recently discussed by other authors [26l|27]. It would be interest¬ 
ing to establish the precise connection between all these constructions. 

Finally, using the notion of fermionic anyon condensation we sketch 
a rough argument demonstrating how one could potentially “simulate” 
a generic fermionic lattice Hamiltonian in 2-|-l dimensions given a copy 
of K^ and a sufficiently rich bosonic system. This argument supports 
the idea that Kd may fully capture the non-local properties of a generic 
fermionic system. 
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1.3 Conclusions and future directions 


Although we have demonstrated the link between fermionic phases 
of matter and spin-TQFTs only in a restricted set of examples, we 
believe that the relation will hold in greater generality. There are two 
natural ways one may try to extend our results 

• It should be possible to adapt our analysis to the general math¬ 
ematical framework of extended topological field theory. The 
analogue of a fermionic phase of matter should be an extended 
topological field theory such that the vector spaces attached to 
{d — l)-dimensional manifolds have Grassmann grading and the 
tensor products are twisted by the Koszul sign rule. The same 
sign combinatorics as in the state sum model should lead to an 
obstruction proportional to [r(; 2 ]- Thus we expect one could prove 
a theorem relating “fermionic” extended topological theories and 
extended spin-TQFTs. 

• It should be possible to give a purely Hamiltonian version of 
our analysis. The fermionic SPT phase associated to the 
z[X^ rj, Pii-i] partition function should give us a recipe to build a 
one-dimensional Hilbert space from a given cocycle Pd-i, which 
could be described as a choice of sign on local patches of the 
space manifold. Hopefully, this recipe will capture the same sign 
ambiguities as one encounters in the construction of a general 
fermionic Hilbert space as a tensor product of local fermionic 
Hilbert spaces associated to local patches of the space manifold. 
Ideally, this would show in full generality, without reference to 
TQFTs, that any fermionic phase of matter can be obtained by 
combining the fermionic SPT phase with some appropriate 
bosonic degrees of freedom and that fermionic phases of mat¬ 
ter should generally require the existence of a spin structure on 
space. 


2 Fermionic SPT phases and the Gu- 
Wen Grassmann integral 


2.1 The Gu-Wen construction 

The standard discrete action for a bosonic SPT phase protected by 
some symmetry group G is built from a U (l)-valued d- cocycle v^{go, • • • , Qd) 



on BG, i.e. a function of {d + 1) G-valued variables, invariant under 
the action of G on itself 

• • • ,ggd) = • • • ,5d) (12) 

and closed under the action of an appropriate differential 

Concretely, given a triangulation of a d-dimensional manifold X 
with a flat G connection, one evaluates the partition function as a 
product over all d-dimensional simplices of evaluated on a lo¬ 
cal trivialization of the connection. The G-symmetry of the cocy¬ 
cle makes the answer independent of the local trivialization and the 
cocycle condition = 1 insures invariance under changes of triangu¬ 
lation. Essentially, the ratio between the partition functions for two 
triangulations which differ by an elementary move equals the partition 
function for the boundary of a (d -|- l)-simplex, which by definition is 
the same as More generally, two triangulations can be related by 
a sequence of moves which can be visualized as a triangulation of a 
cobordism X x [0,1] from the manifold X to itself. 

The Gu-Wen construction of a discrete partition function for fermionic 
SPT phases [7j involves two basic pieces of input: a C/(l)-valued d- 
cochain • • • , gd) on BG and a Z 2 -valued (d—1) cocycle nd-i{go, • ■ ■ , fi'd-i), 
such that 

dl/rf = (13) 

In other words, Vd satisfies the cocycle condition up to signs, which 
are determined from nd -2 through the Steenrod square operation. 

The Gu-Wen partition function can be decomposed into the prod¬ 
uct of three terms, which are not separately invariant under changes of 
triangulation. The first term, which we could denote as or Zy\T\ 
if we want to indicate the specific choice of triangulation T of the 
space-time manifold X, is simply the product over all d-dimensional 
simplices of just as for a bosonic SPT phase. Because Vd is not 
a cocycle, the sign of this term will jump under re-triangulation by 
the Steenrod square of Ud-i integrated over the cobordism from X to 
itself: 

Zy[T] = Zy[T']expm [ Sq^[nd-i] (14) 

dxx[o,i] 

The second term, which we could denote as Zq or Ze\T], contains 
the “fermionic” degrees of freedom. It is a sign, defined by a Grass- 
mann integral whose structure is determined by Ud-i- Schematically, 
one associate a pair of Grassmann odd variables to the two sides of 
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each [d — l)-simplex such that n^-i is 1. The integrand is built as 
a product over d-simplices of the Grassmann variables associated to 
that simplex. We do not expect to be able to write Zg as a standard 
bosonic action, i.e. the integral of the pull-back of some class on BG to 
the manifold: if we could do that, we would have reduced the system 
to a bosonic SPT phase. The Grassmann integral Zg only depends 
on the group variables through the image fdd-i of the cochain Ud-i 
computed on the faces of the triangulation, which is a standard Z 2 co¬ 
cycle. We will see later in Section [6.21 that the Grassmann integral Zg 
coincides with the function <Tn(^, fdd-i) described in the introduction. 

The third term, which we could denote as Zm or Zm[T], is some¬ 
what problematic: it is written in terms of a function md- 2 {go, • • • , gd- 2 ) 
which is not G invariant, but satisfies 6 md -2 = The expression 

for Zm involves a product of (—evaluated over a certain subset 
S of {d — 2)-simphces in T defined in by some local rule: 

Zm[T] = (15) 

ses 

We review the precise definition of S for d = 2,3,4 in section [221 

The product is invariant under re-definitions of md- 2 , but the lack 
of G-invariance of md -2 makes it problematic to define the model 
on a manifold with a non-trivial G-bundle. The wavefunctions built 
from this model have some md -2 dependence which is stripped off 
by hand by the authors of [7] in a non-canonical way, opening the 
possibility for subtle sign changes and ambiguities in the corresponding 
TQFT. Indeed, we will argue here that removing the md -2 dependence 
introduces naturally a dependence on a choice of spin structure on the 
manifold, so that the Gu-Wen fermionic SPT phases are a class of 
invertible spin-TQFTs. 

In order to remove the spurious md -2 dependence, we can imagine 
replacing the collection S of {d — 2)-simphces used in Zm with a col¬ 
lection of [d — 1) simplices E such that dE = S, so that the product 
over S of (—l)™d -2 in Zm can be reorganized to a product Z^ over E 
of 

e£E 

Glearly, this will be only possible if S is exact, and then the result will 
depend on the choice of E: any two different choices of E differ by 
a cycle, which will have in general a non-zero pairing with Ud-i-, and 
thus will give inequivalent partition functions. For a solid torus, this 
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integration by parts should mimick the way the authors of [7] strip off 
the md -2 dependence of wavefunctions. 

At this stage, we have shown that we can construct an improved, 
well-defined partition function if the homology class of S vanishes, and 
that the improved partition function depends on a choice of trivial- 
ization E of S. Although the definition of S depends on the choice 
of triangulation, invariance of the partition function under changes of 
triangulation suggests that the homology class of S should capture 
some intrinsic triangulation-independent property of the underlying 
space-time manifold X. 

We propose that the homology class of S captures precisely the 
second Stiefel-Whitney class [rc 2 ]; which vanishes on spin manifolds. 
We also propose that S itself provides a canonical chain representative 
for W 2 , so that the choice of E actually encodes a choice of a spin 
structure on the manifold. 

The proposal implies that the partition function can only be de¬ 
fined on a spin manifold, and depends on a choice of spin struc¬ 
ture. The combination of signs ZqZ^ will coincide with the function 
z[X, rj , Pd-i] defined in the introduction. Conversely, if the cohomol¬ 
ogy class of W 2 is non-trivial, the partition sum will have an unavoid¬ 
able dependence on re-definitions of md -2 which map to a non-trivial 
{d — 2) cohomology class on the manifold. 

We can offer two strong checks of our proposal: a general consis¬ 
tency check we describe momentarily, and a direct calculation for a 
special class of triangulations, which are obtained by refining a generic 
triangulation by a barycentric subdivision. These triangulations are 
endowed with a canonical representative chain for W 2 , which turns out 
to coincide with 5 for d = 2, 3,4. 

The first consistency check follows from the observation that a 
change of Zm under a change of triangulation is a linear expression 
™'d- 2 [s]) a sum of m evaluated on some collection V oi {d — 2)- 
simplices of the cobordism Y = X x [0,1] from X to itself. The change 
of Zm should only depend on Ud-i, as it is cancelled by the change in 
Ziy and Zq. It should thus be possible to write it as a linear expression 
a sum of Ud-i evaluated on some collection W of (d—1) 
simplices of Y such that dW = V. 

The linear function defined by W must be such that the sum 


^ nd-i[s] + 
sew 



(17) 
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of the variations of and can be cancelled by the variation of 
Ze. 

There is a natural candidate for snch a linear function: the inte¬ 
gral of W2 U Ud-i over T, where W2 is some specific cochain which 
represent the second Stiefel-Whitney class. A neat property of W 2 is 
that Sq‘^[x] -)- rc2 U x is exact for any {d — l)-cocycle x and thus can be 
integrated by parts and plausibly cancelled by the variation of Zq. If 
we write Ud-i = Smd-2 and thus W2 U Ud-i = S{w2 U md-2), we find 
that S must be a representative for W2 on X. 

We can give a direct proof that 5 is a representative for W 2 for a 
special class of triangulations BT, which are a barycentric subdivision 
of some rougher triangulation T. Each d-dimensional simplex s of T is 
subdivided in (d -I- 1)! simplices whose vertices are the barycenters 
of subsets of the vertices of s. More precisely, given a permutation 
a of the vertices of s, the i-th vertex of is the barycenter of the 
vertices (t( 0), • ■ ■ Thus the vertices of have a natural order, 

starting from a vertex of s and ending with the barycenter of s. 

In the Gu-Wen construction, the triangulation is endowed with a 
branching structure, i.e. an orientation of the edges with no closed 
loops, which in turn provides a order to the vertices of simplices: the 
vertex with all outgoing edges is number “0”, the vertex with a single 
incoming edge is number “1”, etcetera. 

If we have a barycentric triangulation, we can simply order the 
1-simplices from the barycenter of fewer vertices to the barycenter of 
more vertices. That gives a useful canonical choice of branching struc¬ 
ture on BT, which induces the same order of the vertices of simplices 
as the natural one of the barycentric subdivision. This branching 
structure has the useful property that each vertex of BT has the same 
position in the order of vertices in all simplices which include that 
vertex: vertices of T are always at position “0” in the order, midpoint 
of segments in T are always at position “1”, etcetera. 

The barycentric subdivision BT has an important property: the 
set of all (d — 2)-simplices provides a chain which is a canonical rep¬ 
resentative for W 2 - The collection 5 of (d — 2) simplices used in [7] 
to define Zm is described as the sum over all (d — 2)-simplices, plus a 
correction term given explicitly in dimensions d = 2, 3,4. We will now 
show by direct inspection that the extra correction term vanishes for 
a barycentric triangulation BT, so that S is precisely the canonical 
representative for 7 x 2 ! 
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Figure 1: The barycentric subdivision of a triangle 

2.2 Barycentric subdivisions in d = 2, 3, 4 

For a general triangulation, the branching structure gives an order to 
the vertices of each triangle: the n-th vertex has n incoming edges. 
The ordering of the vertices also gives an orientation of each simplex. 
This orientation may or not agree with the canonical counterclockwise 
orientation of the simplex. If it does, we have a “+” simplex, if not 
we have a ” simplex. See Figure [2j 

In 2d, the partition function involves a factor of (—for each 
vertex u, and an extra factor of (—foj. each ” triangle which 
has V as the vertex number 1. 

If we pick our triangulation to be a barycentric subdivision BT, 
then each “1” vertex belongs to two — triangles, and we can disregard 
the second contribution. Thus Zm[BT] is the product over all vertices 
V of BT of (—This is simply the pairing of rtii with the 
canonical 0-chain representative S for W 2 - 

Every 2d manifold admits a spin structure, so that W 2 is always 
exact. If we pick some 1-chain E on BT such that dE = S', we can 
replace Zm[BT] with the improved Zj^[BT]. We interpret the choice 
of E as a choice of spin structure on the discretized surface. 

In 2d we can actually sketch a proof that S is a chain represen- 
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Figure 2: The two possible orientations of a triangle and ordering of vertices 
induced by a branching structure on a 2d triangulation 

tative for W 2 even for a generic triangulation T. In 2d, we can build 
representatives for W 2 by taking a vector field and picking the points 
where the vector field vanishes, counting how many times the vector 
field winds around the origin in a neighbourhood of each point, mod¬ 
ulo 2. That is the same as 1 plus the number of times the vector field 
is tangent in the counterclockwise direction to a small circle around 
the point. If we pick the vector field V in Figure [TJ each vertex will 
contribute 1 plus the number of times the vertex appears at position 
“1” in a — triangle. That representative for W 2 coincides with S. 

In 3d, the chain S of edges which appear in Zm consists of all edges 
of the triangulation, together with the (02) edge for all “-t-” tetrahedra 
and the (13) edge for all ” tetrahedra. 

For a barycentric triangulation BT, the (02) edges join vertices of 
T and barycenters of triangles in T. They are shared by two “-I-” trian¬ 
gles. The (13) edges join midpoints of edges in T and the barycenters 
of tetrahedra in T. They are shared by two ” triangles. Thus the 
extra contributions to S cancel out, and S coincides with the sum of 
all edges, i.e. canonical chain representative for W 2 - 

In 4d, the chain S consists of all triangles in the triangulation, 
together with the (013), (134), (123) triangles for a “-I-” 4-simplex, 

(024) for a ” 4-simplex. In order to show these edges are shared 
by an even number of tetrahedra of the corresponding orientation, we 
can observe that they are fixed points of two reflections. 

For example, a (013) triangle has vertices which are barycenters of 
the first, the first 2 and all vertices in a sequence of 4 {vi,V 2 -,v^,V 4 ). 

As it belongs to a face, it is shared by two tetrahedra in T. As is is 
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invariant under exchange of U 3 and U 4 , it is fixed by a second reflection. 

Similarly, a (134) triangle has vertices which are barycenters of the 
first 2, 4, 5 vertices of a sequence of 5 vertices. It is invariant under 
exchanging the first pair, or the second pair. A (123) triangle has 
vertices which are barycenters of the first 2, 3, 4 vertices of a sequence 
of 4 vertices. It belongs to a face and is invariant under exchange of 
the first pair of vertices. Finally, a (024) triangle has vertices which 
are barycenters of the first 1, 3, 5 vertices of a sequence of 5 vertices. 
It is invariant under exchange of the second and third vertices, and 
the fourth and fifth. 

Thus all extra contributions to S cancel out, and S is given by the 
sum of all triangles, i.e. canonical chain representative for W 2 - 

In order to mimic in general dimension d the 2d analysis for a gen¬ 
eral triangulation T, we would need to build a W 2 representative in 
terms of d — 1 vector fields on X. The representative would be concen¬ 
trated at the locus where the vector fields fail to be linearly dependent, 
weighed by the number of times the vanishing linear combination of 
the vector fields winds around the other non-vanishing linear combina¬ 
tions in a neighbourhood of the locus. We expect it should be possible 
to dehne some canonical set of d — 1 vector fields in each simplex, lin¬ 
early dependent at (d — 2 )-simplices, which demonstrate the 5 is a 
representative for W 2 - 


3 The Gu-Wen Grassmann integral as 
a quadratic refinement 

3.1 Quadratic property of the Gu-Wen Grass¬ 
mann integral 

Next, we will look at Zg. We aim to show that the combination 


Z^[BT]Ze[BT] = z[X,r],/3d-i] (18) 

is a quadratic function of the Z 2 cocycle /3d- 1 , which refines the pairing 
on the space of Z 2 (d — l)-dimensional cocycles defined by the higher 
cup product Ud_ 2 . 

The Grassmann integral only depends on the image /3d_i(e) of the 
cochain n^-i on the co-dimension one faces e of the triangulation. 
Thus we can consider some generic (d — l)-cochain /3d_i, assigning a 
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Figure 3: Left: to each edge e such that /9(e) = 1 we assign a 9^ Grassmann 
variable (black dot) and a 6e Grassmann variable (white dot). Each such 
edge contributes dOedOe to the measure. Middle and right: as /3 is a cochain, 
each triangle is associated to two Grassmann variables, which are ordered in 
the integrand according to the grey arrows. 

Z 2 element to each edge of the triangulation, and specialize it to the 
image of Ud-i later on. For notational clarity, we will usually omit 
the subscript {d — 1) from j3. 

In the following formulae we will often refer to a cochain Ap evalu¬ 
ated the simplex defined by vertices oq, • • • simply as vlp(ao • • • Op). 

We will also denote the Grassmann variables associated to an edge 
with vertices a and b as 6ab, Oab- 

We will start with a 2d example, and then proceed to higher di¬ 
mensions. If /3(e) = I for an edge e, we assign to the edge Grassmann 
variables 9e and 9e- 9e is associated to the side where the canonical 
orientation of the face agrees with the orientation of the edge, 9 on 
the side it disagrees. 

To each triangle t we associate an even monomial u{t), the product 
of the Grassmann variables attached to its edges, if present, in the 
order 12, 01, 02 according to the ordering of the vertices for a -|- 
triangle, opposite for a — triangle. 

Thus a -|- triangle has monomial 

./3(12)./3(01)^/3(02) 

'^12 '^01 '^02 

and a — triangle has 

./3(02)^/3(01)^/3(12) 

'^02 '^01 '^12 

We are interested in the sign 

^(/3)= / n d9ed9e n^(t) 

e|/3(e)=l t 


( 20 ) 

( 21 ) 
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We will show directly from the definition that 

a(/3 + /30 = a(/3M/3')(-l)^-^^^' (22) 

i.e. cr(/3) is a quadratic function which refines the intersection pairing 
on 1-cocycles. 

We can combine the two Grassmann integrals in a{l3)a{/3') as 

=/ n dOeSe n d9'je'^'[[u{t)[(3,9,9]uit)[P',e',9'] 

e|/3(e)=l e|/3'(e)=l t 

(23) 

There is a subset of “spurious” Grassmann variables which is associ¬ 
ated to edges for which /3(e) = /3'(e) = 1. If we can systematically 
integrate them out, the remaining Grassmann variables will be as¬ 
sociated to edges for which /3(e) -|- /3'(e) = 1 and coincide with the 
Grassmann variables in the formula for (t(/3 -)- /3'). 

Our strategy is to permute the variables until we form pairs 9e9'f, 
and of spurious variables. These pairs are Grassmann even, 
and can be brought out of the integrand and eliminated against the 
d9ed9ed9'^d9'^ measure to give an overall sign. 

The permutation of variables which brings the spurious pairs to¬ 
gether is very simple: we interleave the variables in the product u{t)\j3, 9, 9]u{t)[f3', 9', 9']: 
for a -|- triangle we have 

= (_ l)/9(01)/3' (12)-H/3(02)/3' (12)-H/3(02)/3' (01). 

• ( 24 ) 

and for a — triangle we have 

= (_ l)/9(01)/3' (02)-H/3(12)/3' (02)-H/3(12)/3' (01). 

• 'eT^ (4)'"'^°'^] (001)'''^°'^] [oT^ (25) 

Thus we have grouped together pairs of spurious variables, if present, 
and left the non-spurious variables in the correct order to be identified 

with u(t)[/3 -)- /3', • ■ ■ 

It is useful to re-distribute the overall sign, by as¬ 

sociating each factor of (—l)^(®)^he) with the triangle to the right of 
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the oriented edge e. Thus each + triangle, which only sits to the right 
of its 02 edge, is associated to an overall factor 

(_1)/9(01)/3'(12)+/3(02)/3'(12)+/3(02)/3'(01)+/3(02)/3'(02) ^ (■_]^)/3(01)/3'(12) 1 ^ 26 ) 

and to each — triangle, which only sits to the right of its 01 and 12 
edges, is associated to an overall factor 

(_1)/9(01)/3'(02)+/3(12)/3'(02)+/3(12)/3'(01)+/3(01)/3'(01)+/3(12)/3'(12) ^ (•_]_)/3(01)/3'(12) 

(27) 

Thus the reorganization of the Grassmann integral a{l3)a{/3') into 
the Grassmann integral cj(/3 + /3') produces a sign of (—l)^(®oi)/3'(ei2) 
for each triangle. This is nothing else but the cup product of f3 and 
/?'! Thus we have the desired quadratic property 

+ = (28) 


We can readily extend the above argument to Zg in any dimen¬ 
sion d. For general d, the Grassmann integral involves again pairs of 
variables associated to the two sides of each oriented (d — l)-simplex 
e such that /3(e) = 1. 

The integrand is a product of monomials u[s] made out of all 
Grassmann variables associated to each d-simplex s. The order of 
the variables in the monomial is determined by a specific rule. We 
will describe the rule later, for now we only need to know that it gives 
a canonical order to the faces of each simplex. 

If we denote as (t(/ 3) again the result of the Grassmann integral, 
the same interleaving operation as in the 2d case gives us immediately 
a proof that cr{/3) is quadratic: 


cr(/?+ /?') = cr(/5)o-(/3')]Je[s,/3,/3'] (29) 

S 


with 


6[s,/ 3,/3'] = (_l)J::re'e./5(‘=)/^'(^')+ESV(e)d(e') 


(30) 


with e > e' in the order determined by u[s] and e > 0 if tt[s, /3] includes 
a 9e variable. We aim to identify the d-cochain in the exponent with 
the higher cup product /3 Ud -2 ■ 

The order of faces induced by u[t] for a -|- simplex places first the 
faces which omit the even vertices (first the one omitting 0, then the 
one omitting 2, etc.), which are labelled by Oe variables, and then the 
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faces which omit the odd vertices (first the one omitting 1, then the 
one omitting 3, etc.), which are labelled by 9e variables. The order for 
a — simplex is the opposite (and the role of Of. and 6 ^). 

We will now compute e[s, /3, /?'] in d = 3 and d = 4 for this choice of 
order and verify it is given by the standard higher cup products /3Ui/3' 
and /3U2/3' respectively. A similar analysis in higher dimension should 
be straightforward. 

Notice that if we were to pick a different choice of order IT for the 
factors in the measure, cTn(/3) would differ from the standard (t(/ 3) 
by a local linear term, which does not affect the quadratic refinement 
property. If we pick a different choice of order IT for the factors in u[t], 
iTn(/3) will differ from the standard (t(/ 3) by a local quadratic term, 
some (/3, /3) pairing defined by a sum over simplices of the product of 
fermion numbers of the variables which have been permuted in u[t] 
to get to the new order. That would change the quadratic refinement 
from /? Urf _2 (3' to some other 

/3Urf_2/3' = /3 U,_2 /3' + (/3, /3') + (/3', (5) (31) 

This is just a different choice of definition for the higher cup product: 
the basic relation in equation 0 remains valid with the re-definition 

/3Urf_3/3' = /3 Urf_2 /?' + 5{(5, (5') + [5(5, (5') + (/3, ,5/3') (32) 

and all other products unchanged. 

3.2 Quadratic refinement in Sd 

At a + simplex, we have an order of faces (123), (013), (023), (012) 
and thus the signs computed from the quadratic pairing 

/3(012)/3'(123) + /3(012)/3'(013) -b /3(012)/3'(023) + /3(023)/3'(123)-b 
+ /3(023)/3'(013) -b /3(013)/3'(123) + /3(023)/3'(023) -b /3(012)/3'(012) = 
/3(023)/3'(012) + /3(013)/3'(123) 

(33) 

At a — simplex, we have an order of faces (012), (023), (013), (123), 
and thus the signs computed from the quadratic pairing 

/3(123)/3'(012) + /3(123)/3'(023) -b /3(123)/3'(013) + /3(013)/3'(012)-b 
+ /3(013)/3'(023) -b /3(023)/3'(012) + /3(013)/3'(013) -b /3(123)/3'(123) = 
/3(023)/3'(012) + /3(013)/3'(123) 

(34) 
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Again, the dependence on the type of triangle drops out. We have a 
d-cochain (/3,/3')3 given by the pairing /3(023)/3'(012) +/3(013)/3'(123). 
This is precisely the definition of /3 Ui /3h Thus we claim that in 3d 
we have 

+ = ( 35 ) 

3.3 Quadratic refinement in Ad 

At a + simplex, we have an order of faces (1234), (0134), (0123),(0234), 
(0124) and thus the signs computed from the quadratic pairing 

/3(0124) [/3'(1234) + /3(0134) + /3(0123) + /3(0234)] + 

+ /3(0234) [/3'(1234) + /3'(0134) + /3'(0123)] + 

+ /3(0123) [/3'(1234) + /3'(0134)] + /3(0134)/3'(1234)+ 

+ /3(0124)/3'(0124) +/3(0234)/3'(0234) = 

/3(0234)/3'(0124) + /3 (0123)/3'(1234)+ 

+ /3(0123)/3'(0134) + /3(0134)^'(1234) (36) 

At a — simplex, we find the same. Thus have a d-cochain 

(/3,/3')4 =/3(0234)/3'(0124) + /3(0123)/3'(1234)+ 

+ /3(0123)/3'(0134) + /3(0134)/3'(1234) (37) 

which is the standard expression for /3 U 2 /3^ 

4 The t’Hooft anomaly 

4.1 Gauge variation of the Gu-Wen Grassmann 
integral 

Next, we can look at the variation of the Grassmann integral under 
exact changes in the cocycle dd-i- The calculation is greatly simplified 
by the quadratic refinement property. 

a(/3+dA) = u(/3)a(<5A)(-l)^v/3Ud-2<5A ^ ^(^)^(5A)(-l)/x/3Ud-3A+Au,_3/3 

(38) 

where we used the basic property of the higher cup product, i.e. eqn. 
m the violation of the Leibniz rule for Ua equals the symmetrization 
of Ua_i. 
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We can specialize the quadratic refinement property to 
fT(5A + 5A') = a((5A)cT((5A')(-l)-^->f 

( 39 ) 

which can be solved up to a linear ambiguity: 

a{5\) = (—1)5Zs 6S ^(^)+/x 

In order to fix the ambiguity, it is useful to compute directly a{5X) 
for the simplest possible case, where A is non-zero only on a single 
{d — 2)-simplex, so that 5X equals 1 on all the {d — l)-simplices whose 
boundary include the selected {d — 2)-simplex. Thus if we go around 
the {d — 2)-simplex, we will encounter a sequence of Grassmann vari¬ 
ables di, with measure factors ±dd 2 idd 2 i+i from the {d— l)-simplices 
and u{t) factors ±d 2 i+ 2 d 2 i+i from the d-simplices inserted around the 
{d — 2)-simplex. 

If all the signs above were -|-, the overall Grassmann integral would 
give a factor of —1. In general, the signs in the measure factors are 
determined by the orientation of the (d— l)-simplices and the signs in 
the integrand factors u[t] are determined by the relative order of the 
two (d — I)-simplices Grassmann variables in u[t]. 

Gonsider for simplicity a barycentric subdivision. We have several 
types of (d — 2) simplices, which can be labellet by two integers from 0 
to d: they are simplices whose vertices do not include the barycenters 
of a vertices or of b vertices. They are associated to an alternating 
sequence of (d—1) simplices which omit either barycenters of a vertices 
or of b vertices. As these simplices are also of alternating type -|- or 
—, it is easy to see that we can pick a direction around the (d — 2)- 
simplex such that the integrand factors will have all -|- signs. Working 
through a few examples, it is easy to convince oneself that the signs 
in the measure factors multiply to 1. Thus if A* is non-zero only on a 
single (d — 2)-simplex s in BT, we can write 

(t((5As) = -1 (41) 

As the quadratic part of equation |40] vanishes for A^, we find that 
S for a barycentric subdivision BT coincides again with the canonical 
representative for W 2 - We expect S to essentially coincide with S and 
the canonical chain representative for W 2 for a general triangulation 
as well. In other words, we expect that 

z[X,r],dXs] = l (42) 
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Thus we can write 


(j((iA) = ^Ud_3<5A+AUd_4A+«)2UA 

and thus: 

_|_ JA) = —l)/x ^'^<i-3'^+'^Ud_3/3+AUd_35A+AUd_4A+-u;2UA 

We can also write the’t Hooft anomaly of z[X, rj, f3d-i] under gauge 
transformations of /3rf_i: 

z[X, r/, (3d-i+6Xd-2] = z[X, q, /3u,_3A+au,_3/3+au,_35a+au,_4A 

(45) 

4.2 A WZW-like expression for a quadratic re¬ 
finement 

For d > 3 one can construct a WZW-like expression for the quadratic 
function of the cocycle /3 as follows. Let us assume that W is a bound¬ 
ary of some compact oriented (d-|- l)-manifold Y. This is automatic if 
d < 3, since the oriented bordism group vanishes for d = 2,3 

m, but in general it is a nontrivial constraint on X. If we are given a 
(d—l)-cocycle /3 G Z'^~^{X, Z 2 ) on X, one can always choose Y so that 
fi extends to a (d — l)-cocycle on Y. To see this, we regard /3 as map 
/3 : X^K{'L 2 , d—1) where K{’L 2 , d—1) is an Eilenberg-MacLane space; 
then the statement we need is that the reduced oriented bordism group 
Q^^{K{’Ii 2 ,d— 1)) vanishes. This follows from the Atiyah-Hirzebruch 
spectral sequence for (unreduced) bordism and the vanishing of the 
reduced homology of K{Z 2 , d — 1) in degree less than d — 1. Moreover, 
if d > 3, the reduced oriented d-dimensional bordism of a product of 
several copies of Ar(Z 2 , d — 1) also vanishes, for the same reason^ This 
implies that one can choose Y to be independent of (3. Such a T is 
not unique, of course. 

Consider now the following WZW-like expression: 

cr(/3) = (-l)/v/5Urf_3/3+«;2U/3 

^This fails for d = 2 , because of the Tor terms in the Kiinneth formula. In fact, 
it is easy to give an example of a closed oriented 2 -manifold X and a pair of classes 
Oi, 02 G H^{X, Z2) such that there is no compact oriented 3 -manifold Y such that X = dY 
and both oi, 02 arise from restriction of cohomology classes on Y. For example, one can 
take A to be a 2 -torus, with oi and 02 being the generators of iJ^(A, Z2). 
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It is independent of the actual choice of Y, as j3 Uds /3 + ^2 U /3 is 
known to be exact if /3 is a cocycle. 

This expression is a quadratic refinement of Dd- 2 - Indeed, if /3 and 
/3' are {d — I)-cocycles on X, we get 


d-{(3 + (3') = d-(/3)5-(/3')(-I)-^x/3Ud_2/3' 


(47) 


It also transforms in the same way as a under 1-form Z 2 gauge 
symmetry: 

a{6X) = (^_l)/y <5^Ud_3<5A-l-«)2U5A ^ 

_ ^'^U,i_45A-l-<5AU,i_4A+AU(^_3(SA-1 -u!2UA 

= (48) 

This means that a and a can only differ by a linear and gauge-invariant 
function of /3. 

Thus for d >3 for all practical purposes we can write 

z[X,r],Pd-i] = (-l)/x^u/3rf_i+Jy/3rf_iUrf_3/3d_i+«>2U/3rf_i ( 49 ) 


5 Fermionic SPT phases and spin cobor- 
dism 

It was proposed in ffH that fermionic Short Range Entangled phases 
in d space-time dimensions with symmetry G and vanishing thermal 
Hall conductivity are classified by the Pontryagin dual of the torsion 
part of Q.^^^{BG). Here BG ~ K{G, 1) is the classifying space of G. 
Some checks of this were performed in [12]. We can now compare with 
the Gu-Wen super cohomology proposal in low dimensions. 

Let us begin with d = 3. In this dimension there are no nontrivial 
fermionic SRE phases in the absence of symmetry, so in the presence 
of symmetry there is no distinction between SRE and SPT phases. 
Prom the mathematical viewpoint, one has = 0, and thus 

Q^pin^BG) coincides with the reduced bordism group Q^^^^{BG). The 
partition function of the model corrected by the spin-structure depen¬ 
dent term is 


Z{X, A, r]) = exp ^ 27 rz j A*z(X, rj, A* 132), 


(50) 
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where M/Z), j32 G Z‘^{BG, Z 2 ), and A is a gauge field on 

X regarded as a map A : X^BG. The cochains /rs and (^2 satisfy 

5z/ 3 = i/?2 U/32. (51) 

It is easy to see that this expression defines an element of the Pon- 
tryagin dual of Indeed, it is clear that Z{A,rj) is multi¬ 

plicative under disjoint union. Now, suppose there exists a compact 
spin 4-manifold Y with boundary (X, ri) such that A extends to a map 
Ay : Y^BG. Then 

exp(27rz / ( 52 ) 

Jx 

On the other hand, the WZW-like expression for z{X,r],l32) becomes 
(_l)/x V'JA*l32+fY W2UA*Y/32+fY A^/32UA^/32 ^ ( 53 ) 

which is clearly the same as (I52p . Since Z(X,A,r]) becomes 1 when 
evaluated on trivial bordism classes, it defines a homomorphism from 
QSpin^BG) to U{1). 

Not all spin cobordism classes can be so obtained. For example, 
for G = Z 2 it is known that Q^^^^{BZ 2 ) — Zg [12], while the Gu-Wen 
construction only gives phases labeled by Z 4 . From the physics side, 
it is also known that 3d fermionic SPT phases with Z 2 symmetry are 
classified by Zg m- 

For d > 3 there may exist nontrivial fermionic SRE phases even in 
the absence of any symmetry. According to m, they exist whenever 
hS'S torsion. If we want to focus on fermionic SPT phases, we 
can restrict to d-manifolds X which define a trivial class in 
Then the same argument shows that if A extends to a compact spin 
{d + I)-manifold Y such that dY = X, then the partition function 
of the Gu-Wen model is I. Therefore, each Gu-Wen supercohomol¬ 
ogy class defines a homomorphism from the reduced bordism group 
Q^^"'{BG) to U{1). Again, in general we do not expect that all such 
homomorphisms can be obtained from the Gu-Wen construction. 

In d = 2 the arguments are a bit different, since there is no 
WZW-like expression for z{X,ri, Pi). The proposal of [H] is that 
fermionic SPT phases are classified by the Pontryagin dual of the 
reduced bordism group Q^^^^iBG). There are also fermionic SRE 
phases in the absence of any symmetry which are classified by the 


24 


dual of = Z 2 . The Gu-Wen construction describes only the 

former. To describe the correspondence, recall that a spin structure 
on an oriented 2-manifold can be identified with a quadratic refine¬ 
ment of the intersection form on Z 2 ) [EldB]. This 

quadratic refinement is nothing but z{X,r], /3i), see Appendix A for a 
detailed discussion. Moreover, it follows from the results of m that 
the value of z{X, rj, /3i) depends only on the bordism class of (A, r], /3i) 
in 0.2^^{ 3 ^ 2 )- Thus the Gu-Wen construction defines a map from 
H^{BG,'Ij 2 ) to spin cobordism of BG. This map is not a homomor¬ 
phism, because z{X,ri,l3i) is not linear but quadratic in /3i. But we 
should remember that every element 1^2 G //^(i?G,M/Z) also gives us 
an element in the spin cobordism of BG. Thus Gu-Wen SPT phases 
are described by pairs ( 2 ^ 2 ,/3i) G H^{BG,M./'Z) x H^{BG,'L 2 ). The 
group structure is a nontrivial extension of H^{BG, Z 2 ) by H^{BG, R/Z): 

(1^2, Pi) + {122, Pi) = (z ^2 + 2^2 + 2^^ P'l^Pi + Pi)- ( 54 ) 

It follows from the Atiyah-Hirzebruch spectral sequence that this ex¬ 
tension is isomorphic to the dual of the reduced bordism ^ 2 ^^{BG), 
in agreement with the proposal of m- 

The nontrivial d = 2 fermionic SRE phase without any symmetry 
is realized by the Kitaev spin chain. The corresponding spin-TQFT 
has the Arf invariant as its partition function (the unique bordism 
invariant of spin structures in d = 2). We recall that the Arf invariant 
of {X,r]) is essentially the average of the z{X,r],Pi) over all Pi G 
Z 2 ) [E]. This spin-TQFT can also be constructed using the 
Gu-Wen Grassmann integral, see section [H 

6 Constructing spin-TQFTs in low di¬ 
mensions 

6.1 State-sum constructions of 2d spin-TQFTs 

An oriented 2d TQFT can be defined axiomatically as a functor from 
a geometric category Cob 2 to the category of vector spaces. The cat¬ 
egory Cob 2 has closed 1-manifolds as objects and oriented bordisms 
between them as morphisms. It is well known that there is a 1-1 
correspondence between oriented 2d TQFTs and commutative Frobe- 
nius algebras. The commutative Frobenius algebra corresponding to 
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a given 2d TQFT encodes 2-point and 3-point functions on a sphere, 
and the rest of the correlators can be reconstructed from it. 

An alternative approach to constructing 2d TQFTs is provided by 
state-sum models [mil]. This approach is more natural from the 
statistical mechanics viewpoint and gives a manifestly local recipe for 
computing the TQFT partition function and correlators for a triangu¬ 
lated closed oriented d-manifold X. One starts with a not necessarily 
commutative semi-simple Frobenius algebra A and defines the parti¬ 
tion function as follows. Fix a basis e*, i G /, in A and denote by 
the structure constants of A in this basis. Let 

9^3 = CkCfi. 

The matrix gij is non-degenerate if A is semi-simple. Let be its 
inverse. It is easy to see that 

Cijk 9il^jk 

is cyclically symmetric. A coloring of a 2-simp lex / of A is an assign¬ 
ment of an element of I to each boundary 1-simplex of /. A coloring 
of a triangulation is a coloring of each 2-simplex. If each 2-simplex 
is colored, each 1-simplex has two colors. The partition function of a 
triangulated manifold A is a sum over all colorings of the triangula¬ 
tion, with the weight of each coloring defined as product of weights 
of 1-simplices and 2-simplices. The weight of a 1-simplex colored by 
i,j G I is 5 *-^. The weight of a 2-simplex whose three edges are col¬ 
ored by i,j,k is Cijk- Here we use the cycling ordering of the edges 
arising from the orientation of A. It is easy to show that the partition 
function thus defined is independent of the choice of a triangulation 

[mnEj. 

The state-sum construction is somewhat redundant, as the par¬ 
tition function and the correlators depend only on the center Z{A), 
which is a semi-simple commutative Frobenius algebra. The algebra A 
can be interpreted as the algebra of boundary operators for a partic¬ 
ular boundary condition for the TQFT based on Z{A). On the other 
hand, the state-sum construction is very explicit and can be easily 
extended to manifolds with boundaries. 

Note that the state-sum construction always gives rise to a semi¬ 
simple Z{A) and therefore does not produce the most general oriented 
2d TQFT. For applications to condensed matter physics, this is not a 
serious drawback, since unitary TQFTs are automatically semi-simple. 
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It is easy to modify the state-sum construction to produce 2d spin- 
TQFTs. One starts with a Z 2 -graded semi-simple algebra A. Let e*, 
i € I, be its basis. Each basis vector is assumed to have a well-defined 
grading /3(i) G Z 2 . As in the bosonic case, we color each edge of each 
2 -simplex with an element of /, so that each 1 -simplex is colored with 
a pair of elements of I. The weights assigned to 1-simplices and 2- 
simplices will be proportional to and Cijk- Since the algebra A 
is Z 2 -graded, the matrix vanishes if /3(z) 7 ^ /3(j)- Thus we may 
assume that if a 1-simplex is colored by z,j, then j3{i) = /3(j). Thus 
each allowed coloring defines a Z 2 -valued 1 -cochain j3. Since Cijk 
vanishes unless /3(z) + (3{j) + (3{k) = 0 , this 1 -cochain is a cocycle. 

In order to account for the Grassmann nature of elements of A, 
we correct the naive sum over the colour of edges by including an 
overall sign. This is the sign which arises from identifying Cijk with 
an element of A* G A* G A* and g"^^ with an element of A ( 8 ) A: we 
pick some order of the factors in each individual weight and then re¬ 
order them to bring together the pairs of A and A* spaces we want to 
contract. 

The combinatorics of the Koszul signs is the same as in the Gu- 
Wen construction. The Grassman variables 9^-, Be act as placeholders 
for the Grassman-odd generators of A* in the weights of 2-simplices 
and the d6e, dBe as placeholders for the Grassmann-odd generators of 
A in the weights of 1-simplices. The contraction between generators 
of A and A* is mimicked by the Grassman integration. The result is 
just cr(/3). 

Thus the naive weight of each allowed coloring will be the product 
of g^^ over 1-simplices, Cijk over 2-simplices and cr{f3). We can rewrite 
the sum over colorings as a sum over colorings producing a particular 
1-cocycle (3 followed by a sum over /3. Let us denote by Z[/3] the result 
of the first summation. The discussion in section 2 implies that Z[/3] is 
independent of the triangulation up to a sign. To make it completely 
independent, we need to choose a trivialization rj of W 2 and multiply 
Z[/3] by a correction factor 

i.e. we take the correct weight to be the product of g^^ over 1 -simplices, 
Cijk over 2-simplices and z{X, g, (3). Then the partition partition func¬ 
tion depends on the spin structure 77 on A, but not on a particular tri¬ 
angulation. Notice that in 2d z{X, g, (3) coincides with the well-known 
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quadratic refinement of the intersection form, evaluated on 

(3 (see Appendix A). 

The simplest example of a 2d spin-TQFT is obtained if we take 
A to be the Clifford algebra Cl(l). It is generated by 1 and an odd 
variable rj satisfying 7^ = 1. The 2x2 matrix gij, i £ Z 2 , is given by 
Qij = 26ij, while Cijk is equal to either 2 or 0 depending on whether 
the sum of the indices is 0 or 1 modulo 2. In this case the coloring is 
completely determined by the 1-cochain /3, and the weight is simply 
z{X,r],/3) = (—More precisely, if we denote by E and F 
the number of 1-simplices and 2-simplices respectively, the partition 
function reduces to 

2 F-E = 2i-5WArf(A,r/), 

P 

where g{X) is the genus of A, and Arf(A, g) E {±1} is the Arf invari¬ 
ant of the spin structure g [EKinilSS]: 

Avf{X,g) = 2 -sW E 

[p]&H^{X,Z2) 

The partition function of the state-sum is not ±1, but it differs from 
it by an exponential of an integral of a local counter-term (the Euler 
density), and thus can be made ±1 by a local redefinition. After such 
a redefinition, we get an invertible 2d spin-TQFT which describes the 
basic fermionic SRE phase in two space-time dimensions (the May- 
orana spin chain). This spin-TQFT has also been briefly discussed in 

m- 


6.2 State-sum construction of 3 d spin-TQFTs 

In 3d, the state-sum construction of bosonic TQFTs is known as the 
Turaev-Viro construction [T9l |^. Again, it does not produce the 
most general oriented 3d TQFTjj, but it has the advantage that it can 
be easily extended to manifolds with boundaries. 

The starting point of the Turaev-Viro construction is a spherical 
fusion category A. It can be thought of as a categorification of a finite¬ 
dimensional semi-simple algebra A. A fusion category is a semi-simple 
rigid monoidal category with a finite number of simple objects. Let 

®For example, Chern-Simons TQFTs for simple Lie groups do not arise from the Turaev- 
Viro construction. 
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I be the set of simple objects. The partition function is a sum over 
colorings of l-simplices of the triangulation with elements of I. If we 
denote such a coloring by </>, we can write 

Z = '^w[4>\Z^ (55) 

(t> 


The weights w[(l)] are a product over 1-simplces of an appropriate 
function of their color, times a product over all O-simplices of an ap¬ 
propriate ^-dependent constant. 

Each individual term in the partition function can be thought 
of as the evaluation of a tensor network: each face with edges of color 
i,j,k is associated to a pair of vector spaces Vijk, each one asso¬ 
ciated to a side of the face, and each tetrahedron t to chosen element 
in the tensor product of the four vector spaces associated to its 
faces. The tensors are contracted together at the common faces. The 
partition function will be invariant under changes of triangulations if 
the tensors satisfy some basic axioms, such as the pentagon 

axiom, which guarantees invariance under a 2 — 3 move. 

The same construction applied to a manifold with a boundary 
produces wave-functions, which can be reproduced as ground states 
of a Hamiltonian built from the same data m- 

The fermionic version of this construction is given in the literature 
in the “Hamiltonian” form, without a description of the correspond¬ 
ing partition partition function. In this section we will attempt to 
assemble the ingredients of the fermionic construction into a partition 
function. 

The starting point is a spherical super-fusion category A. As far 
as we know, this notion was first mentioned in [9] and appeared more 
recently in [25] in the context of fermion condensation (see below). 
We will now describe what this means in a somewhat informal way, 
relegating the details to Appendix C. A super-fusion category has a 
finite number of simple objects Vj, i G I representing independent 
quasi-particle excitations. There are no nonzero morphisms between 
Vi and Vj if i 7 ^ j, while the space of morphisms from Vj to Vj is a 
Z 2 -graded division algebra (i.e. every nonzero element has an inverse). 
Unlike in the bosonic case, there are two possibilities for such a division 
algebra: C or a Clifford algebra with one generator Cl(l). Following 
[9| we will call the former kind of simple object a bosonic simple 
object, and the latter kind of object a Major ana simple object. Since 
Cl(l) ( 8 ) Cl(l) ~ Cl( 2 ), the identity object is always bosonic. 
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This can be explained in a different way. One can tensor any object 
of A with a Z 2 -graded vector space. Tensoring with a bosonic vector 
space of dimension N is equivalent to taking N copies of the quasi¬ 
particle. Tensoring with (a fermionic vector space of dimension 1) 
has the physical meaning of attaching a fermion to the quasi-particle. 
The latter operation maps an object to an isomorphic one, but the 
obvious isomorphism is odd (fermionic). There may or may not be an 
even isomorphism as well. Now, suppose V) is a simple object. If V) 
is even-isomorphic to ( 8 ) Vj, then Hom(Vj,V)) ~ Cl(l), i.e. Vi is 

Majorana. On the other hand, if Vj is not even-isomorphic to C^I^( 8 )Vi, 
then Hom(yj,V)) ~ C, i.e. Vi is bosonic. 

As usual, one can fuse quasi-particles, and the fusion rules are 
described by the way tensor products of simple objects decompose 
into sums of simple objects; 

Vi ( 8 ) Vj ~ 8 ) Vfc- (56) 

Here the “coefficients” H^j are finite-dimensional Z 2 -graded vector 
spaces. If some of the simple objects are Majorana, there is a subtlety: 
the absolute grading on some of is not well-defined. For simplicity, 
we will assume that all simple objects are bosonic, so that all have 
a well-defined grading. 

The fusion of quasi-particles is associative, in the sense that for 
any i,j, k ^ I there is an even isomorphism 

a.{i,j, k) : (V) 8 Vj) 8 Vk^Vi 8 {Vj 8 14). 

Expanding both sides in terms of simple objects, we deduce that 
a{i,j,k) is determined by an even linear map of Z 2 -graded vector 
spaces 

pU j I 

k m n 

This map is known as a 6 j symbol. The collection of all 6 j symbols 
satisfy an associativity constraint which ensures that the two ways of 
constructing an isomorphism from (((V) 81^)814) 8 V;) to (Ej8(V)8 
(I48V;))) are the same. To write down a concrete form for it, we need 
to choose a basis , aV g in each space Hfp Let e(aV) be the 
fermionic parity of the vector e^k . Then the associativity constraint 
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(a.k.a. the fermionic pentagon equation) reads: 
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( 57 ) 


Note the sign on the r.h.s. It reflects the fact that composition of mor- 
phisms in the symmetric tensor category of Z 2 -graded vector spaces 
is dehned using the Koszul sign rule. 

A spherical structure on a super-fusion category assigns to every 
object V its dual V*, so that V** ~ V and that for every two objects 
U, V we have isomorphisms 


Hom(C/, V) ~ Hom([/ ® F*, 1) ~ Hom(y* 0 [/, 1) ~ 

~Hom(l,C/*(8)y) ~Hom(l,y®[/*). (58) 

Let us denote by the dual of a simple object V). The spherical 
structure ensures that the spaces Hijk = Hfj are cyclically symmetric 
and that the space is dual to Hijk for all j, k € I. 

In the fermionic Turaev-Viro construction one considers colorings 
of 1-simplices of an oriented triangulated 3-manifold X with objects 
Vi- We can choose a branching structure on the triangulation, so that 
the vertices of each 3-simp lex and each 2-simp lex have an order. To 
each 2-simplex one can attach a Z 2 -graded vector spaces Hijk and its 
dual, depending on the orientation. Next we would like to attach a 
6j symbol, which is an element of a tensor product of four Z 2 -graded 
vector spaces, to a 3-simplex, and form the partition function by con¬ 
tracting these tensors along their shared 2-simphces. But there are 
two problems with this. First, the 6j symbol itself is an element of the 
tensor product with a particular ordering of Z 2 -graded factors H^j and 
their duals. Changing this order may change the sign of some com¬ 
ponents of the 6j symbol, but there is no natural order on the faces 
of a 3-simplex. Second, the contraction of two odd elements of dual 
Z 2 -graded vector spaces changes sign if one exchanges their order, but 
for a given 2-simplex in an oriented triangulation there is no preferred 
order of two 3-simplices sharing it, unless the 2-simplex itself is given 
an orientation. Such an orientation is induced by a branching struc¬ 
ture, but this means that the partition function might depend on the 
branching structure. 
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These choices are analogous to the choices made in the construction 
of fermionic SPT phases: the order of the factors in the tensor prod¬ 
uct defining 6 j symbols is analogous to the order of the Grassmann 
variables in u[t] and the choice of order of the two sides of a 2 -simplex 
is analogous to the choice of dOdO or dOdO in the fermionic measure. 
One can make the analogy completely precise in the following way. 
First, we choose a basis 


^ Jijkj 

in every space Ffjjfc.The weight of a particular coloring can be com¬ 
puted by summing over the variables aijk attached to every 2 -simplex 
whose edges are colored by i,j,k. If we assume as before that the 
basis vectors have definite parity, each choice of a variables gives a 
Z 2 -valued 2-cochain. Further, since 6 j symbols are even maps, this 
2 -cochain becomes a 2 -cocycle (32 when restricted to configurations of 
a variables for which all 6 j symbols are nonzero. If we perform the 
summation over such variables first by fixing (32, and then summing 
over all choices of (32, then the Koszul sign is simply ctii{(32), where 
n denotes the ordering choices made in the definition of the partition 
function. 

Thus the tensor network amplitude restricted to some fermionic 
parity sector (32 can be decomposed as 

Z4/3) = Zj(/32)an(/32) (59) 

where is a “bosonized” amplitude computed without keeping track 
of Koszul signs in manipulating the vector spaces in the tensors, but 
rather pretends that all vector spaces are bosonic. 

There is also a close analogy between the sign factor (—l)^(“U"i)e( 4 i 9 ) 
in the fermionic pentagon equation and the twisted cocycle condition 
satisfied by the 3-cochain used in the construction of 3d fermionic 
SPT phases. Indeed, we can write the twisting factor as@ 

= ('_i)/ 32(012)/32(234) ^ ('_^y/32U/32](01234) ^-gQ^ 

®In order to write this formula, we have assigned a specific order 01234 to the vertices 
of 4-simplex, associated respectively to edges imnp, ijs, kjmq,knl, Ipqs. This order is 
the one associated with the branching structure induced by the orientation of edges in 
the definition of the 6j symbols, such that a factor of is associated to a positively 
oriented face, with edges i and j positively oriented and k negatively oriented, and {H^)* 
is associated to a negatively oriented face, with edges i and j negatively oriented and k 
positively oriented 
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This is precisely the cup square of ^2 evaluated on the 4-simplex asso¬ 
ciated to the 2-3 movelU Notice that the fermionic pentagon relation 
involves a sum of terms with different fermionic parity in the inter¬ 
nal faces, but the sign factor only involves the fermionic parity of the 
external faces. 

The fermionic pentagon relation is written in terms of the com¬ 
ponents of the tensors, and thus does not include any Koszul signs. 

The sign factor describes a change in the sign of the bosonic weight 
under a 2-3 move. If we put the Koszul signs back in, the change of 
<711 (/32) under a 2-3 move will almost cancel the sign coming from the 
pentagon equation: it will trade the cup square of /32 with the usual 
linear term f W2(52 controlled by the canonical representative for W 2 - 
Again, this linear term will provide an obstruction to assembling a 
well-defined partition function, unless the underlying manifold admits 
a spin structure. 

It should be now clear how to assemble a well-defined spin-structure- 
dependent partition sum: 

??, h] (61) 

where we decomposed the partition sum over eigenspaces of /32, de¬ 
noted the “bosonized” partition sum in each sector as Z^[f32] and 
combined the Koszul signs and the spin-structure-dependent correc¬ 
tion into z[X,r],j32]- The change of z[X,r],j32] under 2-3 moves will 
cancel against the sign factor in the pentagon relation, giving a well- 
defined object. 

Notice that Z[r]] is almost written as the partition function of a 
topological Z 2 gauge theory with a 2-form gauge field (32, except that 
Z^[/32] is not well-defined on its own, as the 2 — 3 move involves a 
summation over both bosonic and fermionic vector spaces when gluing 
the tetrahedra. Inspection of the 2 — 3 move, though shows that 
the cocycle condition for (32 almost fixes the Grassmann parity on 
internal faces, up to a binary choice, corresponding to a shift of (32 by 
the coboundary of a cochain concentrated on the the “t” edge of the 
tetrahedra on the 3 side of the 2 — 3 move. B 

^We would like to thank Zhengcheng Gu for pointing it out. 

®Notice that because of that, the correction J^rjU P may add an extra relative sign 
between terms of a 2 — 3 move, essentially controlled by the cup product between W 2 and a 
cochain concentrated on the “t” edge of the 4-symplex. That relative sign does not affect 
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That means that we can define a sensible bosonic partition function 
by summing only over /?2 which differ by an exact cocycle 

z‘[&i = EE + ai](-l)/v/32UAi+AiU/3.+AiU5Ai 

Ai 4 > 

(62) 

and re-write the fermionic partition function as 

Z[r?] = ^Z'[/32 ]z[X,7?,/32] (63) 

02 

Thus the fermionic Turaev-Viro theories can be obtained by gaug¬ 
ing the diagonal 1-form Z 2 global symmetry of a product of appro¬ 
priate bosonic theories and the spin-TFT defined by the partition 
function z[X,ri,(32]- 

7 Fermionic anyon condensation 

7.1 Self-fermions and anomalous l-form Z 2 global 
symmetry 

In this section we would like to argue that in 2 -|- 1 dimensions, a 
bosonic TQFT equipped with a 1-form Z 2 global symmetry with the 
same’t Hooft anomaly as is the same as a TQFT which has in the 
spectrum of quasi-particles a self-fermion e which fuses with itself to 
the identity. 

Our first observation is that a TQFT which has in the spectrum 
of quasi-particles a self-boson b which fuses with itself to the iden¬ 
tity is automatically equipped with a non-anomalous 1-form Z 2 global 
symmetry, whose generators are simply b worldlines. The operation 
of gauging the 1-form symmetry is the same as the operation of anyon 
condensation applied to b. 

If we are given any two TQFTs with self-fermions e and e' respec¬ 
tively which fuse with itself to the identity, the product of the two 
theories has a self-boson ee'. This shows that if we take / and f 
worldlines to be generators of 1-form Z 2 global symmetries, the two 
Z 2 global symmetries will have the same Z 2 -valued’t Hooft anomaly. 

the 2 — 3 move with the standard branching structure, where the “t” edge has vertices 
(13) and does not contribute to a cup product. We expect that relative sign to be present 
in the pentagon relation for other choices of branching structure. 
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That means we only need to compute the ’t Hooft anomaly for 
some representative TQFT with a self fermion e which fuses with 
itself to the identity. All other TQFTs with such a quasi-particle will 
have the same’t Hooft anomaly. 

A simple choice of such a TQFT is the toric code, aka Z 2 gauge 
theory. We can describe it by the discrete partition function 

(_l)/x“iU'56i (04) 

with oi and bi being discrete cochains. 

The e,m and e quasi-particles correspond to Wilson loops for oi, bi 
and ai-|- 6 i respectively. The former two quasi-p^articles are generators 
of two non-anomalous 1-form Z 2 symmetries. □ For example, we can 
couple the theory to a Z 2 2-cocycle /3| as 

(_l)/x'UU<5fei-HaiU/3| (00) 


The action is invariant under gauge transformations /3| -|-(5Af, accom¬ 
panied by bi bi + Xf. Similarly, we can couple the theory to a Z 2 
2 -cocycle /3™ as 

(_l)/x“iu<5bi+/3rufei ( 00 ) 

The action is invariant under gauge transformations -|- bXf', ac¬ 
companied by ai — 7 > oi -|- AT*. 

On the other hand, the 1-form Z 2 symmetry associated to e is 
encoded in the action 


(_l)/x 


aiU(5fei+aiU/32+/32Ufei 


(67) 


Under gauge transformations (52 + (5Ai, accompanied by 61 ^ 61 -|- Ai 
and ai — 7 > oi -|- Ai, the action varies by 




AiU(5Ai+AiU/32+/32UAi 


( 68 ) 


This is precisely the anomalous variation for K^\ This verifies our 
claim. 


7.2 Examples of fermionic anyon condensation 

We will come back to the toric code momentarily. Before that, we 
should give a simple examples of bosonic and fermionic theories related 
by fermionic anyon condensation. 

®The two symmetries have a mixed anomaly. 
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Consider first a U{l)ik Chern-Simons theory. The Wilson loop of 
charge 2k fuses with itself to the identity. If k is even it is a self-boson 
b, while if k is odd it is a self-fermion e. For even k, condensing b leads 
to a U{l)k CS theory. We expect the same to occur for odd k, but 
then the resulting U{l)k CS theory is a spin-TQFT. This is consistent 
with the fact that this is a fermionic anyon condensation. 

It is interesting to observe that the toric code has a Z 2 global 
symmetry which exchange the e and m particles. This symmetry is 
rather hard to make explicit in a concrete description of the model. 
For example, the action we wrote above is not invariant under a Z 2 
symmetry transformation exchanging ai and bi: 


(-l)i'v 


aiUSbi +6iU(5ai 


-^\Jx 5aiUi5fei 


(69) 


Assuming that the Z 2 global symmetry is not broken by the fermionic 
condensation of e, the resulting spin-TQFT should still enjoy the Z 2 
global symmetry. Indeed, there are good reasons to believe this is 
a fermionic SPT phase, the root fermionic SPT phase which is not 
captured by the super group cohomology classification. For example, 
we can consider the Hilbert space of the fermionic theory on a general 
Riemann surface (equipped with spin structure). Sectors with a non¬ 
trivial 2-form gauge field flux are described by a Riemann surface with 
a single e insertion, which do not contribute ground states. It is quite 
reasonable to assume that the 2^ ground states with no /32 background 
flux connection will be projected down to 2^ one-dimensional Hilbert 
spaces labelled by the choice of spin structure. 

A Z 2 -invariant boundary condition for the toric code has gapless 
modes; a non-chiral 2d Ising model coupled to the toric code in such 
a way that e ends on the spin operator a{z,z), m on the dual spin 
operator fi{z, z) and e ends on chiral or anti-chiral fermion operators 
i^{z) and ^(z), while the energy operator e{z,z) can exist without 
a quasi-particle world line ending on it. The Z 2 global symmetry is 
Kramer-Wannier duality. 

After condensing the e quasi-particle, the chiral or anti-chiral fermion 
operators ipiz) and V’(^) in the boundary theory will be able to be 
inserted in correlation functions without a quasi-particle world line 
ending on it. Thus the boundary theory reduces to a theory of free 
Majorana fermions, one chiral and the other anti-chiral, with the Z 2 
global symmetry acting only on one of the two fermions. This is an 
the expected feature of the root fermionic SPT phase with Z 2 global 
symmetry. 
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Another interesting example is the fermionic spin-TQFT “s-Ising” 
associated to the Ising 3d TQFT. The Ising TQFT has three quasi¬ 
particles; 1, a and e. In the presence of a boundary, these quasi¬ 
particles end on the corresponding operators of a chiral Ising model. 
We expect that upon condensing the e quasi-particle, the correspond¬ 
ing chiral operator on the boundary, which is a free fermion operator, 
will be free to appear by itself. Thus the boundary theory should 
consist of a single chiral Majorana fermion, which can indeed be well- 
defined on a manifold equipped with a spin structure. 

It is also interesting to look at the Hilbert space of s-Ising on 
Riemann surfaces of various genus, equipped with a spin structure. On 
a torus, the Ising TQFT has three states. Upon gauging the 1-form 
symmetry, we expect them to be projected to three one-dimensional 
Hilbert spaces labelled by the three even spin structures on the torus. 
As they have no /32 hux on the Riemann surface, they should have 
even fermion number. Upon gauging the 1-form symmetry we also 
gain twisted sectors with a non-zero flux of /32 on the Riemann surface 
and odd fermion number, i.e. a single insertion of e. On the torus, 
that gives an extra state (labelled by a loop of a with an e insertion), 
which we expect to give a one-dimensional Hilbert spaces labelled by 
the single odd spin structure on the torus, of odd fermion number. 

On a genus 2 Riemann surface the Ising TQFT has 10 untwisted 
sectors, which should map to 10 one-dimensional Hilbert spaces la¬ 
belled by the 10 even spin structures on the torus. We can also find 6 
states for the Ising TQFT on a genus two surface with an extra e inser¬ 
tion, which should map to 6 one-dimensional Hilbert spaces labelled 
by the 6 odd spin structures on the torus. 

In general, we expect the s-Ising spin-TQFT to have one-dimensional 
spaces of ground states for every choice of spin structure on a Riemann 
surface, with even fermion number if the spin structure is even and 
odd if the spin structure is odd. In other words, if we compactly the s- 
Ising TQFT on a circle with Ramond boundary conditions, we should 
be left with a 2d theory which computes the Arf invariant of Riemann 
surfaces. 

It is somewhat surprising to observe that the s-Ising spin-TQFT 
has one-dimensional spaces of ground states, and yet supports chiral 
edge modes and would likely not be considered a fermionic SPT phase. 
Indeed, the standard classihcation of fermionic SPT phases in 2 -|- 
1 dimensions with fermion number symmetry only predicts no non¬ 
trivial theories. 
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7.3 Bosonization in 2 + 1 dimensions 

The operation of fermionic condensation on the toric code is somewhat 
intriguing from the perspective of defining a notion of bosonization in 
2 + 1 dimensions. Suppose that we were given a version of the toric 
code lattice Hamiltonian with the property that creating fermionic 
quasi-particle excitations could be created at much lower energy cost 
than the e and m quasi-particles. The Hilbert space of the low-energy 
theory will be the even part of a fermionic Fock space, but the inter¬ 
actions will not be completely local: fermion bilinear operators would 
still have to be represented by open string operators, corresponding 
to an e worldline connecting the two fermionic creation/destruction 
operators. 

Suppose also that we could tensor that system with it's and gauge 
the diagonal 1-form Z 2 symmetry. Then the e worldline string oper¬ 
ators, dressed by the self-fermion quasiparticles of ifs, would become 
invisible in the gauged theory. This would liberate the endpoints of 
the open worldline string operators, which would behave as fermionic 
creation or destruction operators. The Hilbert space of the gauged 
low-energy theory should become a fermionic Fock space, which one 
could use to simulate some general fermionic Hamiltonian in a local 
way. 

It would be interesting to make this recipe more precise. Ideally, 
one should find a lattice Hamiltonian description of the toric code 
which both allows the fermion quasiparticles to proliferate and which 
can be coupled appropriately to a /?2 cocycle. 
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A The Gu-Wen Grassmann integral in 
2d as a quadratic refinement 

Notice that we can depict /? as a collection of closed, disjointed paths 
Ca in the dual graph to the triangulation, sequences of edges e“ with 
l3{e^) = 1 sharing a common triangle with the previous and next 
edges in the path. We can order the monomials in sequences along 
each path, and factorize the Grassmann integral as 


o-(/5)= n 


]ldGetd0etllu{ei] 


(70) 


where we ordered the integrand in one direction and the measure in 
the opposite direction. We assumed that e[t,l3] is 1 unless a triangle 
is crossed by some path Ca- 

For each path Ca, denote as c“ the sequence of Grassmann variables 
encountered along the path. We have 

/ JJc2j_|_iC2i+2 = —1 (71) 

i i 


where ia is the length of Ca- 

The integration variables in d6e°-d6e9 are encountered almost in 
the order dc^ij^idc^i, up to permutations of two variables associated 
to the same edge. Similarly, the integration variables are encountered 
in the product WiU{C_i) almost in the same order as Hi'^ 2 i+i'^ 2 i+ 2 , 
up to permutations of two variables associated to the same triangle. 

Thus the overall sign of the integral can be computed by multiply¬ 
ing a sign for each edge and each triangle along Ca'- 


o-(/5)= n 

a 




(72) 


where the edge sign ei is 1 if the edge orientation points to the right 
of Ca, —1 if it points to the left of Ca, while the triangle sign 62 is 1 if 
both edges of the triangle crossed by Ca point to the right of Ca, — 1 
if both point to the left, while if they point in opposite directions the 
sign is 1 if Ca is turning left at the triangle, —1 if turning right. 

It is useful to combine the sign associated with crossing an edge 
to the sign associated to transversing the subsequent triangle. The 
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Figure 4: A collection of edges with /9(ei) = 1, organized into a path in the 
dual graph to the triangulation. The Grassmann variables are encountered 
along the path in an order which differs by simple local permutations from 
the order they appear with in the measure and integrand of the Grassmann 
integral. 

+ 

0 0 




Figure 5: As a path C proceeds along the triangulation, it accumulates a 
sign when crossing each edge, and a sign when passing through each triangle, 
depending on if the Grassmann variables it encounters are set in canonical 
order in the integral, or not. In this figure we indicate the directions along 
which a path receives a + sign. The opposite directions produce a — sign. 
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Figure 6: As a path C proceeds along the triangulation from just before 
crossing an edge to just before crossing the next, it accumulates a sign as 
indicated in this hgure: +1 in alms all cases, except the ones indicated by a 
dashed line. 

+ 

0 0 




Figure 7: A useful reference vector field V inside a triangle equipped with 
branching structure. 

resulting sign is almost always +1, except when we cross a 01 edge 
and turn towards a 12 edge. See Figures [5] and [H 

There is a simple interpretation of such sign. Consider a certain 
canonical vector field V inside each triangle, nonzero away from the 
vertices and continuous across the edges, which flows out of the vertex 
0 and into the vertex 2, turning around near the vertex 1. This can 
be thought of as a continuation inside the triangle of the branching 
structure directions, as in Figure [71 

Given a smooth, non-self-intersecting path C in the triangulation, 
which does not pass through the vertices, we can measure how many 
times the vector held V restricted to C winds around compared to the 
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tangent vector of C. We can measure such winding number w{V,C), 
modulo 2, by counting how many times the vector field V becomes 
tangent to the path C. Inspection of the sign rule of[6l shows that the 
overall sign accumulated by C is (— 

Thus we conclude that 

a 

if /3 is represented by the collection of non-intersecting paths Ca- 

The sign (t{(3) does depend on the choice of cochain /3 in a given 
cohomology class: if we sweep some path C across a vertex v of the 
triangulation, the winding number wiV, C) will generically jump. It 
is natural to expect, and it should be easy to show, that the jump 
modulo 2 will be the same as the contribution of that vertex v to the 
set S of vertices which contribute to Zm, i-e. 1 plus the number of 
times V is the vertex 1 of a — triangle. 

Indeed, that would insure that 

a(/3 + 5A)=a(/3)n(-l)^^"^ (74) 

v£S 

SO that the variation of Zq = cr{n) cancels against the variation of Z^ 
under mi —>■ mi + A. 

This fact is clearly true if the triangulation is a barycentric sub¬ 
division BT: the 0 vertices are sources, the 1 vertices simple saddles 
and the 2 vertices are sinks. All these configurations are associated to 
a jump of the winding number by one unit. 

The combination 

JJ (_i)/3(e) (75) 

e£E 

is a known expression for the quadratic refinement of the intersection 
form on Z 2 1-cocycles on a surface [23]. Besides transforming in the 
obvious way under changes of spin structure E, it satisfies 

QeW + /3') = QEiP) + QEiP') + 113 up' (76) 

Thus 

Z^[BT]Zg[BT] = (77) 

If we have a generic triangulation, we can still try to assemble a 
smooth vector field V and look at the zeroes of V, with appropriate 
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multiplicity, as a representative chain S for W 2 - In general, we expect 
solutions of dE = S to still describe a choice of spin structure on 
the discretized surface. One may interprete measuring the winding 
numbers w(V, C) modulo 2, corrected by the pairing with E, as a way 
to build a spin bundle on the triangulated surface. 

B Higher cup products 

As described in the introduction, higher cup products satisfy the re¬ 
cursive property 

AUa B + B A = 5{A Ua+i B) + 6A Ua+i B + A Ua+i 6B (78) 

with Uq = U. 

There is a rather explicit canonical formula for the higher cup 
products [23]: 

[Ap Ua Bq] (0, • • • , P + (?) = 

^ - , 004 , ••• 079 ) 




where the sequences of arguments of Ap and Bq must agree with the 
degree p,q of Ap and Bq. 

For example, we recover the standard cup product: 


[ApU Bq]{0,--- ,p + q) = A{0,---p)B{p,--- ,p + q) (80) 


The next cup products are 


[ApUiBq]{0,--- ,p+q-l) = '^A{0,---io,io+q,p+q-l)B{io,--- ,io+q) 


(81) 


and 


[Ap U 2 Bq] (0, • • • ,p + q — 2) 



jo<n 


• 5(io, • • • ,n,p-l- ii - io - 1, • • • ,p + q-2) (82) 


Etcetera. 
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C Super-fusion categories 

A category C enriched over Vectk (the category of finite-dimensional 
vector spaces over a field k) is a category whose morphism sets Hom(A, i?) 
are vector spaces over k, and the composition of morphisms is bilinear. 

It can be thought of as a generalization of the notion of an algebra 
over a field (“algebra with many objects”). Similarly, a category C 
enriched over sVectk (the category of Z 2 -graded vector spaces over a 
field k) is a category whose morphism sets Hom(A, i?) are Z 2 -graded 
vector spaces over k, and the composition of morphisms is bilinear and 
is compatible with the grading. Such a catregory C can be thought of 
a generalization of the notion of a Z 2 -graded algebra over a field. In 
what follows we will suppress k (in physical applications one typically 
has k = C). 

An important difference between ordinary algebras and Z 2 -graded 
algebras is the way the tensor product of algebras is defined. For 
ordinary algebras, the product of algebras C and D is an algebra 
whose underlying set is C 0 D and the multiplication is defined by 

(c 0 d) • (c' (g) d') = cc' (g> dd', c, c' G C, d, d' € D. 

On the other hand, for Z 2 -graded algebras one defines the multiplica¬ 
tion by 

(c (g) d) • (c'(g> d') = (—• c') (g> (d • d'), c,c'&C, d,d' € D, 

where e(a) G Z 2 is the fermionic parity of a G C. This is the Koszul 
sign rule. This has an analog for categories: the Cartesian product of 
categories over sVect is defined differently from the Cartesian product 
of categories over Vect. In the former case, if C, C', C" and D, D', D" 
are objects of C and T> respectively, the composition of morphisms in 
C X T> is defined as follows: 

(a 0 /3) • (7 (g) d) = (—(a • 7 ) 0 (/3 • d), 

a G Homc(C',C "),7 G Homc(C, C'), 

13 G Homp(T)',T»")>'5 G Rom^iD, D'). (83) 

while in the latter case the sign factor on the r.h.s. is absent. 0 

^°More generally, one can consider categories enriched over a symmetric tensor category. 
The definition of Cartesian product for such categories depends on the symmetric tensor 
category. 
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This in turn has consequences for how a tensor product on a cate¬ 
gory is defined. Both in the ordinary and the Z 2 -graded case, a tensor 
product on C is a functor m from C x C to C, m : {X,Y) (->■ X (g) T, 
VX, Y € Ob(C), together with an associator. The associator is a natu¬ 
ral isomorphism between two functors from CxCxC toC built from m. 
The first one is —), —) and the second one is m(—, m(—, —)). 

The associator must satisfy the pentagon equation. Since the defini¬ 
tion of the Cartesian product depends on whether we consider cate¬ 
gories over Vect or sVect, so does the definition of the tensor structure. 

Let us specialize to the case of (multi)-fusion categories. Multi¬ 
fusion categories are rigid (and in particular semi-simple) monoidal 
categories (categories with a tensor product and a unit object) over 
Vect or sVect which have finitely many isomorphism classes of simple 
objects. To emphasize the difference between the two cases, we will 
refer to multi-fusion categories over sVect as super-multi-fusion cat¬ 
egories, while multi-fusion categories over Vect will be simply called 
multi-fusion categories. While multi-fusion categories are module cat¬ 
egories over Vect, super-multi-fusion categories are module categories 
over sVect. 

Let C be a super-multi-fusion category. For every object X G 
Ob(C) we have a Z 2 -graded algebra Homc(V, V). If X is simple, this 
algebra must be a division algebra, therefore it is isomorphic either to 
C or Cl(l). In the former case we will say that X is a bosonic simple 
object, while in the latter case we will say that X is a Majorana simple 
object. 

A super-fusion category is a super-multi-fusion category whose unit 
object 1 is simple. Since 1 (g) 1 ~ 1, this implies that Hom(l, 1) = C, 
i.e. 1 is a bosonic object. 

Since a super-fusion category is rigid, for every object V we have 
a dual object V*, and even morphisms evL ■ V* (g) V—7>1, coev£, : 
l^V 0 V*, evij : V 0 V*^l and coevjj : l^V* 0 V satisfying the 
usual identities. One can show [21] that V** is isomorphic to V. 
A pivotal structure on a super-fusion category is a choice of such 
isomorphisms for all V in a way compatible with the tensor product. 
Given a pivotal structure, one can define the left and right dimensions 
of V by composing coev^ and evR or cobvr and ev^. The left and 
right dimensions are complex numbers. If the left and right dimensions 
are equal for all V, one says that the pivotal structure is spherical. 
A spherical super-fusion category is a super-fusion category equipped 
with a spherical pivotal structure. 


45 


References 


[1] X. Chen, Z. -C. Gu, and X. -G. Wen, “Local unitary transforma¬ 
tions, long-range quantum entanglement, wave function renor¬ 
malization, and topological order,” Phys. Rev. B 82, 155138 
( 2010 ). 

[2] A. Kitaev, talk at the IPAM workshop “Symmetry and Topology 
in Quantum Matter,” Jan 26-30, 2015. 

[3] E. Witten, “Supersymmetric index of three-dimensional gauge 
theory,” In *Shifman, M.A. (ed.): The many faces of the super¬ 
world* 156-184 |hep-th/9903005|. 

[4] X. G. Wen, “Quantum field theory of many-body systems: Prom 
the origin of sound to an origin of light and electrons,” Oxford, 
UK: Univ. Pr. (2004). 

[5] X. Chen, Z. -C. Gu, and X. -G. Wen, “Glassification of gapped 
symmetric phases in ID spin systems,” Phys. Rev. 83, 035107 
( 2011 ). 

[6] L. Fidkowski, A. Kitaev,‘ ‘Topological phases of fermions in one 
dimension,” Phys. Rev. B 83, 075103 (2011). 

[7] Z-C. Gu and X-G. Wen, “Symmetry-protected topological or¬ 
ders for interacting fermions: Fermionic topological nonlinear a- 
models and a special group supercohomology theory,” Phys. Rev. 
B 90, 115141 (2014) [arXiv:1201.2648 [cond-mat]]. 

[8] Z-C. Gu, Z. Wang, X-G. Wen, “Lattice model for fermionic toric 
code,” Phys. Rev. B 90, 085140 (2014) [arXiv: 1309.7032 [cond- 
mat]]. 

[9] Z-C. Gu, Z. Wang, X-G. Wen, “A classification of 2D fermionic 
and bosonic topological orders,” Phys. Rev. B 91, 125149 (2015) 
[arXiv:1010.1517 [cond-mat]]. 

[10] M. A. Levin and X. G. Wen, “String net condensation: A Physi¬ 
cal mechanism for topological phases,” Phys. Rev. B 71, 045110 
(2005) |cond-mat/0404617|. 

[11] A. Kapustin, “Symmetry Protected Topological Phases, 
Anomalies, and Cobordisms: Beyond Group Cohomology,” 
arXiv: 1403.1467 [cond-mat.str-el]. 

[12] A. Kapustin, R. Thorngren, A. Turzillo and Z. Wang, “Fermionic 
Symmetry Protected Topological Phases and Cobordisms,” 
arXiv: 1406.7329 [cond-mat.str-el]. 


46 


[13] Z-C. Gu and M. Levin, “The effect of interactions on 2D 
fermionic symmetry-protected topological phases with Z 2 sym¬ 
metry,” Phys. Rev. B 89, 201113 (2014) [arXiv: 1304.4569 [cond- 
mat]]. 

[14] R. E. Stong, “Notes on cobordism theory,” Princeton University 
Press (1968). 

[15] M. F. Atiyah, “Riemann surfaces and spin structures,” Ann. Sci. 
Ecole Norm. Sup. 4, 47-62 (1971). 

[16] D. Johnson, “Spin structures and quadratic forms on surfaces,” 
J. London Math. Soc. 22, 365-373 (1980). 

[17] M. Fukuma, S. Hosono and H. Kawai, “Lattice topological field 
theory in two-dimensions,” Commun. Math. Phys. 161, 157 
(1994) [hep-th/9212154|. 

[18] C. Bachas and P. M. S. Petropoulos, “Topological models on the 
lattice and a remark on string theory cloning,” Commun. Math. 
Phys. 152, 191 (1993) [hep-th/9205031 . 

[19] V. Turaev and O. Ya. Viro, “State sum invariants of 3-manifolds 
and quantum 6j symbols,” Topology 31 (1992), no. 4, 865902. 

[20] J. Barrett and B. Westbury, “Invariants of piecewise-linear 
3-manifolds,” Trans. Amer. Math. Soc. 348 (1996), no. 10, 
39974022. 

[21] V. Drinfeld, S. Gelaki, D. Nikshych, V. Ostrik, “On braided fu¬ 
sion categories,” Selecta Mathematica 16, no. 1, 1-119. 

[22] G. W. Moore and G. Segal, “D-branes and K-theory in 2D topo¬ 
logical field theory,” hep-th/0609042. 

[23] D.Cimasoni and N. Reshetikhin, “Dimers on Surface Graphs and 
Spin Structures,”. Comm. Math. Phys. 275 (2007), no. 1, 187-208 
math-ph /0608070. 

[24] N. Steenrod, Products of cocycles and extensions of mappings, 
Ann. of Math. (2) 48, 290 (1947). 

[25] K. Walker, talk at the IPAM workshop “Symmetry and Topology 
in Quantum Matter,” Jan. 26-30, 2015. 

[26] J. W. Barrett and S. O. G. Tavares, “Two-dimensional state sum 
models and spin structures,” arXiv:1312.7561 [math-ph]. 

[27] S. Novak and 1. Runkel, “State-sum constructions of two- 
dimensional topological quantum field theories on spin surfaces,” 
arXiv: 1402.2839 [math.QA]. 


47 




